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We show using numerical simulations that a rich variety of novel colloidal crystalline states are 
realized on square and triangular two dimensional periodic substrates which can be experimentally 
created using crossed laser arrays. When there are more colloids than potential substrate minima, 
multiple colloids are trapped at each substrate minima and act as a single particle with a rotational 
degree of freedom, giving rise to a new type of orientational order. We call these states colloidal 
molecular crystals. A two-step melting can also occur in which individual colloidal molecules initially 
rotate, destroying the overall orientational order, followed by the onset of inter-well colloidal hopping. 
PACS numbers; 73.50.-h 



Colloidal particles are an ideal system for studying 
2D ordering and melting as the individual particle posi- 
tions and dynamics can be directly visualized. Colloidal 
melting has been extensively studied for systems with 
smooth substrates, where evidence for defect-mediated 
melting transitions has been found Another av- 

enue of study has been the examination of colloidal crys- 
tallization and melting in the presence of a ID periodic 
substrate, which are typically created with laser arrays 
iQ-^. These studies have revealed a novel laser- induced 
freezing, as well as a reentrant laser induced melting as 
a function of substrate strength . A far less studied 
case is colloidal crystallization on 2D periodic substrates. 
Here numerous commensuration effects can be expected 
when the periodicity of the colloidal lattice matches with 
the substrate lattice constant. Several recent experimen- 
tal studies have demonstrated that it is possible to create 
a 2D substrate for colloids using optical tweezer arrays 
I pot , templating and 2D crossed laser arrays [ p^ . 

Almost nothing is known, however, about what type of 
colloidal crystalline states could occur and what types 
of melting can occur in this system. Although colloids 
interacting with ID substrates have been investigated 
theoretically [^,0 and numerically the case of a 2D 
substrate has not. The study of colloidal crystals with 
2D substrates may also lend insight into other condensed 
matter systems which can be described as elastic parti- 
cles on 2D substrates. Examples of such systems include 
atoms on atomic surfaces [p^ and vortices in supercon- 
ductors with periodic pinning arrays [ p^ . It should also 
be possible to realize entirely new states in the colloidal 
system, since when there are more colloids than substrate 
minima, each minima will contain multiple colloids which 
can have their own internal degrees of freedom. In addi- 
tion, the control of different types of colloidal crystal ge- 
ometries can be of technological importance for the con- 
struction of certain devices, such as photonic band-gap 
materials. 

In this work we examine colloidal ordering and melt- 
ing on 2D periodic square and triangular substrates using 
Langevin simulations. We find that a rich variety of novel 



colloidal crystalline states can be achieved. We focus on 
the case of filling fractions at which there are an integer 
number of colloids per substrate minima. When there are 
more colloids than potential substrate minima, multiple 
colloids are trapped at each minima. The colloids within 
each of these minima can act as a single particle with in- 
ternal degrees of rotational freedom, forming trimer and 
dimer states that have an additional long-range orienta- 
tional order with respect to neighboring minima. These 
states are similar to molecular crystals; hence, we term 
them colloidal molecular crystals or CMC's. These crys- 
talline states also exhibit a novel multi-stage melting 
where, for low T, the colloidal molecules possess both ori- 
entational and translational order. At higher T, the ori- 
entational order is lost as the colloidal molecules in each 
minima began to rotate, but individual colloidal diffusion 
does not occur and the colloids are still confined in each 
substrate minima. For large enough temperature, indi- 
vidual colloidal diffusion occurs. We map out the phase 
diagram for temperature vs substrate strength and find 
that two stage melting only occurs for sufficiently strong 
substrates, and that the transition temperature from the 
solid to the disordered solid in the presence of a sub- 
strate is lower than the same transition temperature for 
a sample without a substrate. In addition, for fixed tem- 
perature and increasing substrate strength, we observe 
a transition from the ordered solid phase to a partially 
ordered solid phase, similar to the reentrant melting seen 
for colloids interacting with ID periodic substrates. 

We simulate a 2D system of Nc colloids with peri- 
odic boundary conditions in the x and y directions, us- 
ing Langevin dynamics as employed in previous colloidal 
simulations |§,|l5| . The equation of motion for a colloid i 
is: 

f^f« + f, + fr (1) 

Here fy = — X^^j ViV^(rij) is the interaction force from 
the other colloids. The colloid-colloid interaction be- 
tween colloids i and j is a Yukawa or screened Coulomb 
potential, V{rij) — (Q-^/jri — |) exp(— K|ri — rj|), where 
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FIG. 1. The colloid configurations (black dots) at T = 0.0 
for a square 2D periodic substrate with A = 2.5, for differ- 
ent densities of colloids with a fixed substrate period, (a) 
One colloid per potential minima with the colloids forming a 
square commensurate lattice, (b) Two colloids per minima 
form a colloidal dimer state, with each dimer perpendicular 
to neighboring dimers. (c) Three colloids per minima form a 
trimer state with orientational ordering, (d) Four colloids per 
minima produce an aligned quadrimer state. 
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FIG. 2. The colloidal configuration (black dots) at 
T = 0.0 after annealing a triangular 2D periodic substrate 
with A = 2.5 for an increasing density of colloids and a 
fixed substrate period, (a) One colloid per substrate minima 
forms a commensurate triangular lattice, (b) Two colloids 
per minima form an orientationally ordered dimer state, with 
the same dimer orientation in every other row. (c) Three 
colloids per minima form an orientationally ordered trimer 
state. (d)Four colloids per minima form a pattern of aligned 
diamonds. 



Q is the charge of the particles, 1/k is the screening 
length, and Fj^) is the position of particle i (j). The 
length of the system is measured in units of the lattice 
constant ao and wc take the screening length 1/k = ao/2. 
For the force from the 2D substrate, fs, we consider both 
square and triangular substrates with strength A, pe- 
riod fixed at flo, and Nm minima. For square substrates, 
fg = Asin(27ra;/ao)x + Asin(27rj//ao)y, and for triangular 
substrates, fs = Y^i=i ^sin(27r]3j/ao)[cos(0i)x— sin(0j)y], 
where pi — xcos{Oi) — t/sin(^?;) + ao/2, 9i = 7r/6, 
02 = 7r/2, and 6*3 = 57r/6. The thermal force fr is a ran- 
domly fluctuating force from random kicks. We start the 
system at a temperature where all the colloids are diffus- 
ing rapidly and gradually cool to T = 0.0. We define the 
melting temperature at zero substrate strength as T^. In 
this model, we do not take into account hydrodynamic ef- 
fects or possible long-range attractions between colloids. 
We have conducted a wide variety of simulations for dif- 
ferent system sizes, ao values, and screening lengths 1/k, 
and find that all qualitative results are robust. 

In Fig. 1 we show the colloidal positions for a system 
with a square substrate. In Fig. 1(a) the colloidal period- 
icity matches one to one with the substrate periodicity, 
= ^m, so that each colloid is located at the center of 



the potential minima and a square colloidal crystal forms. 

In Fig. 1(b), for twice as many colloids Nc = 2Nm, each 
minima now captures two colloids. Neither colloid is lo- 
cated at the exact minima; instead, they are ofii'set from 
it due to the colloid-colloid repulsion. The states inside 
the minima can be regarded as a colloidal dimer with a 
rotational degree of freedom. Fig. 1(b) also shows that 
there are f,wo types of ordering occurring: one arising 
from the square substrate, and the other due to the spe- 
cific rotational orientation of the colloidal dimers, with 
neighboring dimers perpendicular to one another. The 
orientational ordering of the dimers is due to the colloidal 
repulsion, and allows the distance between the colloids 
to be maximized imder the constraint of the square sub- 
strate. If the colloid interaction range is short, such that 
one dimer pair does not interact with the neighboring 
pair, the orientational dimer ordering seen in Fig. 1(b) 
does not occur. In Fig. 1(c), the colloidal configuration 
for Nc = 3Nm again shows two types of ordering, with 
each minima capturing three colloids that form a trimer 
state. The trimers have a specific orientational order in 
which trimers in adjacent columns of minima are rotated 
by 60°, due to the repulsion of neighboring trimers. In 
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FIG. 3. The colloidal positions (black dots) and trajecto- 
ries (lines) over fixed time intervals at different temperatures 
for the case of the dimer state on a square substrate shown 
in Fig. 1(b). (a) T/T^ = 0.25, where is the melting tem- 
perature at zero substrate strength. The system is in the 
ordered solid phase and only small vibrations of the colloids 
around their equilibrium positions occur, (b) T/T^ = 1.5. 
The dimers rotate within the substrate minima, destroying 
the orientational order of the dimers, but the colloids remain 
trapped inside each minima so the system is still in a solid 
phase, (c) T/T^ = 4.0. Individual colloidal diffusion occurs 
throughout the sample and the system is in the liquid phase. 

Fig. 1(d), A^c = 4:Nm also has an ordered superlat- 
ticc structure with orientational ordering. The colloidal 
quadrimer states are all aligned in the same direction. 
For increasing colloidal densities Nc = nNm one can 
expect even more ordcrings with each potential minima 
capturing n colloids. Since the colloids in each minima 
act as a single object with an internal degree of freedom, 
similar to a molecular dimer or trimcr, and since these 
states are confined in a lattice, we call these states col- 
loidal molecular crystals (CMC's). 

In Fig. 2 we illustrate the CMC's on a triangular sub- 
strate. For one colloid per minima, Nc = Nm, shown 
in Fig. 2(a), the colloids form a commensurate triangu- 
lar lattice with each colloid located at the center of the 
potential minima. For = 2Nm, a dimer state occurs 
in each minima as seen in Fig. 2(b). As in the square 
pinning case, an additional orientational ordering of the 
dimers occurs. The dimers in each row have the same 
orientation, which is rotated 45° with respect to the ad- 
jacent rows. For Nc = 3Nm, illustrated in Fig. 2(c), a 
timer state forms and an orientational ordering occurs in 
which, unlike the square case [Fig. 1(c)], all the trimers 



have the same orientation. In Fig. 2(d) a superlattice 
state also occurs with the colloid molecules forming a 
diamond substructure with a triangular superstructure. 

We find that it is also possible to create partially or- 
dered and disordered states at fractional filling fractions, 
such as 1.5 colloids per minima, Nc = 1.5Nm. In this 
case, for the square siibstratcs every other minima forms 
a colloidal dimer state; however, the dimers do not have 
long range orientational order. Similarly, for the trian- 
gular substrates it is not possible to arrange a state such 
that minima with only one colloid will be between every 
other dimer, so again dimer ordering is not observed at 
fractional filling. 

In Fig. 3 we show a representative example of the two 
stage melting of the CMC, in the dimer state for the 
square substrate [Fig. 1(b)]. As illustrated in Fig. 3(a) 
at low temperatures T/T^ < 0.25 (where is the melt- 
ing temperature at zero substrate strength), both orien- 
tational and translational order are present and the sys- 
tem is frozen. We label this phase the ordered solid. In 
Fig. 3(b), at T/T^ = 1.5, the dimers begin to rotate 
within the minima; however, diffusion of individual col- 
loids throughout the sample does not occur so the system 
is still frozen with the dimer orientational order lost. We 
label this phase the partially ordered solid. For a higher 
temperature T/T^ — 4.0 [Fig. 3(c)], the system enters 
a modulated liquid phase. Here the colloids began to 
diffuse throughout the system. We have observed this 
same type of multi-stage melting for the other CMC's in 
Figs. 1 and 2. 

In Fig. 4 we show the phase diagram of the tempera- 
ture T/T^ vs substrate strength A for the case of a square 
substrate and Nc — 2iVm. The melting line is determined 
from the onset of diffusion, and the solid to disordered 
solid transition line is determined from the correlation 
between the dimers as well as the onset of dimer rota- 
tion, measured using the average particle displacements. 
For zero siibstrate A = 0, the colloids form a triangular 
lattice and the clean melting temperature is deter- 
mined from the onset of particle diffusion. For increas- 
ing A the melting temperature monotonically increases. 
The ordered solid melts directly to the liquid for A < 2.0. 
For A > 2.0 an intermediate rotational melting transition 
between the ordered solid and partially ordered solid oc- 
curs. The rotational melting temperature is always less 
than the clean melting line, since both these melting tem- 
peratures are determined by the elastic properties of the 
colloids, which is maximum for a triangular lattice in the 
zero substrate limit. The ordered solid is not triangular 
and hence should have less elasticity. In addition, the 
phase diagram shows that for increasing A the transi- 
tion line from the ordered solid to partially ordered solid 
decreases in temperature. This occurs since, as A is in- 
creased, the dimer length decreases, the dimers become 
effectively further apart, and the dimer-dimer interaction 
strength goes down. This result also implies that at fixed 
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FIG. 4. Phase diagram of temperature T/T^ vs substrate 
strength, for the square substrate with two colloids per poten- 
tial minima. is the melting temperature for a sample with 
no substrate. Black circles: The melting line as determined 
from the onset of diffusion. Open circles: The transition to 
the partially ordered or rotating dimer phase as determined 
from the orientational correlation of the dimers. 

temperature it is possible to observe a transition from 
the ordered to partially ordered solid by increasing the 
strength of the substrate. This is similar to the reen- 
trant melting observed for increasing substrate strength 
for colloids in ID periodic substrates j^Q. We find quali- 
tatively similar phase diagrams for square and triangular 
substrates for both the dimer and trimer states. For the 
trimers the onset of rotational melting occurs at weaker 
substrate strengths than for the dimers, due to the de- 
creased anisotropy of the trimer state compared to the 
dimers. In addition for stronger screening the transition 
from the ordered solid to the partially ordered solid oc- 
curs for both lower temperature and substrate strength. 
An interesting issue is the nature of the melting transi- 
tions, which is beyond the scope of this work. 

To summarize, we have demonstrated that a rich va- 
riety of novel colloidal crystalline states can be achieved 
with square and triangular two dimensional substrates. 
We find that when there are more colloids than substrate 
potential minima, certain colloidal molecular crystals ap- 
pear in which each minima captures an equal number of 
colloids. The colloids in these traps can act as a single 
particle with a rotational degree of freedom, such as a 
dimer or trimer. We also show that CMC states exhibit 
a multi-stage melting for sufficiently strong substrates, 
where the orientational order of the colloidal molecule 
states is lost first followed by the translational order. We 
map out the phase diagram and show that the solid to 
disordered solid transition occurs for temperatures lower 
than the clean melting temperature. It is also possible to 
have a transition at fixed temperature for increasing sub- 
strate strength from a ordered solid to a partially ordered 



solid, similar to the reentrant melting observed for col- 
loids on ID periodic substrates. Since the colloids within 
a minima can act as a single particle with a rotational 
degree of freedom, our results also suggest that certain 
canonical statistical mechanics models, such as Ising, XY, 
Potts, and frustrated models, may be realized with col- 
loids on two dimensional periodic substrates. The states 
predicted here should be observable for colloids interact- 
ing with crossed laser arrays or optical tweezer arrays, 
dusty plasmas in 2D with periodic potentials, and vor- 
tices in superconductors with periodic substrates. 
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